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The stochastic quantisation technique of Parisi and Wu is extended to study non-equilibrium 
statistical mechanics. We show that this scheme is capable of handling white as well as coloured 
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A quarter of a century back , Parisi and Wu [1] pro- 
posed a method of quantisation based on a stochastic 
Langevin dynamics of a physical system in a fifth time 
r. They showed that at the perturbative level, the usual 
quantum field thery (at the level of the correlation func- 
tions) was recovered in the r — > oo of this dynamics. 
In this work, we extend this technique to study non- 
equilibrium statistical mechanics. 

To study systems out of equilibrium, the usual practice 
is to use response fields [2]. These are used to determine 
propagators and correlators. As these are unrelated, one 
has to separarately calculate them. 

In the stochastic quantisation scheme propagators and 
correlation functions in the fictitous "time" dimension 
are related in the usual way. Moreover, the equilibrium 
limit , viz. t — > oo gives back the real time propagators 
and correlation functions. This way of looking at the 
problem allows us the luxury of a fluctuation-dissipation 
in the fictitous time and what is more important, our 
procedure works unaltered even when the real time noise 
source is coloured. 

1. Parisi- Wu Stochastic Quantisation 

We begin with a brief review. Euclidean quantum field 
theoretic correlation functions for the field 6 correspond- 
ing to an action S[<p] are given by [1,3] 

< Q\T6(x 1 )6(x 2 )...6(x l )\0 > 



^ D6[6{x 1 )6{x 2 )..6{x l )] e - S ^ 



(1) 



Parisi and Wu proposed the following alternative 
method: 

(a) Introduce an extra fictitous "time" dimension r 
in addition to the four spacetime x M and postulate a 
Langevin dynamics of the field 6 in the "time" r as 



d6(x,t,r) SS 



<9r 
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+ r](x,t,r) 



(2) 



where r\ is a gaussian random variable satisfying < 
r/(x,t,T) > v = 0; < r){x,t,T)r)(x' ,i' ,t') > n = 25 (x — 
x')S(t-t')S(T-T') 



(b)Evaluate the stochastic average of the fields 6 V sat- 
isfying (2) viz. 

< 6{xi,Tl)6{x 2 ,T 2 )...6{x h Tl) > v 

(c)Put T\ = t 2 = ...Ti and take the limit t\ — > oo. Then 
one has 

lim Tl ^oo < <j>(Xl,Tl)<f)(Xl,Tl)...<l>(xi,Tl) >,, 



/ D</»[</)(x 1 )0(x 2 ).. ( /»(x i )]e- s ^ 



jDtj* 



-sr. 



(3) 



Thus the real time correlation functions are obtained in 
the equilibrium limit of the fictitous time r.Why this is 
possible is understood from the associated Fokker-Planck 
dynamics. The probability of finding the system in the 
configuration 6 at time t is P(6, t) and the evolution of 
P in time t is governed by the Fokker-Planck equation 



dP 
~dt 



S 2 P 



86 \ 56 J 



(4) 



The above equation can be recast into Schroedinger-type 
equation f| = -2H FP ty with f = P(6, t)e s ^/ 2 and 



H FP = -(1/2)- 



8S\ 2 



MS) 



r> 2 <? 



(5) 



The Fokker-Planck hamiltonian H FP is positive semi- 
definite with ground state energy zero and ground state 
*o = e" 5 ^]/ 2 and Um t ^aoP(<l>,t) = constant e ~ s[4,]/2 . 
2. The Edwards- Wilkinson (EW) Model 
One of the earliest models for surface growth is the lin- 
ear Edwards- Wilkinson model [4] described by the equa- 
tion of motion 



dh(x, t) 
dt 



vV 2 h(x,t) +/3(x,t) + c 



(6a) 



where h is the surface's position,/? the stochastic white 
noise and c the constant average velocity of the propa- 
gating surface./? satisfies 

< /3(x 1 ,t 1 )j3(x 2 ,t 2 ) >= 2eS(x 2 -x 1 )S(t 2 -t 1 ) (66) 
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It can be shown that corresponding to the Langevin equa- 
tion 



dh(x,t) SS 



dt 



Sh 



+ (3(x,t) 



(7) 



the Hamiltonian Hew = / dx(Vh) 2 gives back equa- 
tion (6a). We apply the technique first to the Edwards- 
Wilkinson model which is a linear theory.The model is 
described in real time t by the equations (7) or (6). 
The probability distribution associated with (3 is 



P{(3) OC e ~i;fd D xdtf3( X ,t) 2 = e -±fd D xdt(!£-vV 2 h) 2 



(8a) 



= e 



where 



s = hf (0j^(- 2 + ^ 4 )Mfc,-)M-fc,--) (86) 

Introduce the fictitous time "r" such that 
dh{k,io,r) SS 



dr 



Sh 



+ T](k,LO,' 



(9a) 



< r](k, u, r)r/(fc', J, r') = 25(t - r')5{k + k')S{uj + uj') 

(%) 

(96) ensures that the "fluctuation-dissipation" theorem 
is satisfied in r(cvcn if it may not hold in real time t) 
and that as r — > oo, P(k,uj,r) — > P(/3) ,the equilibrium 
distribution. To calculate the propagator, there are two 
approaches: 

(a) From r-variable Langevin Equation 

The equation is 

dh(x,t,r) _ 2 



vV l h{x,t,T)+ri{x,t,T) (10) 



The solution in momentum space is 

f-OO 

h{k,r)= dT'e- vk2{ - T - T '^{k,T') 
Jo 

and the propagator is 

< h(k,u>,T)h(—k,—(j,T') 

= "T2 [e vk2]r ~ T ' ] - e-^ 2 l T+r 'l]<5(r - r' 



(11) 



In the limit of r = r' oo all the usual EW results are 
obtained. 

(b) From the Fokker-Planck action 

In the functional integral approach to stochastic quan- 
tisation [3] , one has a " euclidean" lagrangian L FP corre- 
sponding to H FP 



dh 



L" = (l/2)l = l +lm 



l?V-n/4)fS 



\5h 2 



(12) 



and the Fokker-Planck action is 

S FP = J dx dt drL FP (13a) 

In (k,uj,uj T ) space, this action takes the form 
oFP ' _ [ ,dk du) du) r 

h{k, lu, lu t ) {^lo 2 t + +[(lu 2 + v 2 ^)] 2 ^ h(-k, -uj, lo t ) (136) 

Thus the height-height correlation function in configura- 
tion space is 



< 



n,,„ / i \ f°° , dk . .duo . ,duj T . 

h(x,t,T)h(0,t',T') >= J j — )(—)(^) 



■ikx 



<j(t—t ) — iLO T {r — r') 



LO 2 + {iO 2 + t/ 2 fc 4 ) 2 



(14) 



In our treatment the Langevin dynamics is effectively 
taken to start at r = — oo so that at any finite time 
t, which is infinitely removed from the starting time, 
we are already in the equilibrium distribution. Putting 
t = t' ,t = t' and completing the integrations give 



< h(0)h(x) >-x 1 



(15) 



so that the index a — 1/2. 

We now consider the scaling behaviour for the height 
variable h. Assume ,for some length scale L, that t ~ 
L z , h ~ L a in D space dimensions. Note that we regard t 
as a parameter only and hence its scaling properties are 
meaningless. This is more so justified as t — > oo gives 
us the usual real time correlation functions. It readily 
follows that a = 1 — D/2, z = 2 for L FP as in equation 
(17). So for D = 1 , a = 1/2. 

Note that if we had started from a modified Langevin 
equation of the form. 



dh\ = SS 
dt \Sh 



+ v 



the only change would have been a factor P in front of 
the term proportional to ( ^ttt I in the action. But this 




term does not contribute to the model considered here 
and so there will not be any change in our conclusions. 

3. Edwards- Wilkinson model with with 
coloured Gaussian noise 

For coloured Gaussian noise (in real time) the EW 
model is described by the equations 



dh(x, t) 
dt 



= vV 2 h{x,t) + (}{x,t) 



(16) 
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with ft satisfying (17). The action is 




(—itj + ^fc 2 )p.(k + p)/i(p, u')h(— k — p, — to — uj')h{k,uj) 



FIG. 1: 
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< j8(fci,wi)/9(fe 2 ,W2) >= + fc 2 )%i+w 2 ) (17) 



Writing the fictitous time Langevin equation as 

where 77 satisfies (96). The Fokker-Planck equation corre- 
sponding to (9a) is 



d 2 P 



d -l = ^—(P dS ) + 

<9r dh(k,u) dh(—k,—w) dh(—k,w)dh(—k,u)) 

The Fokker-Planck hamiltonian is now 

So by going over to the fictitous time , it is always possible 
to write a Fokker-Planck hamiltonian even if the original 
equation in real time does not have the structure |£ = 

Proceeding as before the propagator in coordinate 
space in the limit n = t 2 — > 00 is 



< 



f dke lkx f dw dw T 



LU 2 + [k 2 P(0J 2 + v 2 k 4 )] 2 



A+ 2 p 



Therefore the index a = 1/2+ p. Now consider the model 
of Medina et al[ref.6] 

9h{x ^ T) = uV 2 h(x,t,r) {\/2){Vhf +P(x, t) (21) 



-Xk 2p 



d p d q du)\ du>2 
(2w) D (2tt) d 2tt 2tt 



[p.(k - p)][q.(k - q)]/i(p,wi)/i(k - p,w- Wi) 

ft(q, w 2 )/i(k - q, w - w 2 )^ (22) 

In fictitous time the Langevin eqution is 
dh{k,u),r) SS 



dr 



Sh(—k, —lj, t) 



+ r?(fc,w,r) (23a) 



< r?(fc, u, T) V (k', u', t') = 25{t - r')S(k + k')5(uj + uj 1 ) 

(236) 

(236) ensures that the "fluctuation-dissipation" theorem 
is satisfied in r.The evolution to the equilibrium limit is 
governed by S FP constructed from the action S as in 
(22) and using (12) and (13). The fluctuation-dissipation 
theorem relates the connected two point function to the 
imaginary part of the response function i.e. 



C = —ImG 



(24) 



where C refers to connected part of the two-point func- 
tion in fictitous time r and G refers to the full response 
function. The full response function has the structure 



G 



-i0J T + k 2 P(uJ 2 + k 2z ) 



while 



C = 



k 4 P(u> 2 + k 2z ) 2 )+u 2 



Upto one loop level the relevant diagrams are in fig.l. 
Of these the first one is like a mass insertion vertex (of 
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perturbation order O(A)) and hence ignorable. However, 
the second one (0(A 2 )) depends on the external mo- 
menta. Now considering upto one loop level 

G-i „ G -i(0) _ s 

where £ contains the relevant one loop contribution, i.e. 
the second diagram. Denoting by V the contribution 
from the vertex 



d D p dui duj T 
(2tt) d (2tt) (2tt) 



VVGC 



Now consider the momentum behaviour. Noting that 
vk 2 ~ k z , the second diagram gives S behaviour as ~ 



G- 1 ~ k 2 P+ 2z .The vertex V" behaves as - k 2 Pvk 2 .k 2 ~ 
fc 2p+z+2 , while dimensional matching implies that u> ~ 
fc z . Therefore fc 2p+2z = k D+z+4 P +2z+4 - i P- 4z . This gives 
z= i(£> + 4-2p). 

This result agrees exactly with what is found in refs. 
(4) and (5). Since we expect results to be more singular 
with p > 0, we expect the validity of our calculation 
only for that range of p for which z is smaller than the 
corresponding z for the KPZ problem. For D = 1, this 
restricts the validity of our analysis to p > 1/4, which 
is exactly what was found in ref.(5). In summary, then 
we have a technique where the out-of-equilibrium system 
with a coloured noise can be handled by techniques of 
equilibrium statistical physics. 
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